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Abstract 

The phase diagram of the two-dimensional extended one-band U-V-J Hub- 
bard model is considered within a mean-field approximation and two- and 
many-patch renormalization group (RG) approaches near the van Hove band 
fillings. At small t' and J > mean-field and many-patch RG approaches 
give similar results for the leading spin-density-wave (SDW) instability, while 
the two-patch RG approach, which predicts a wide region of charge-flux (CF) 
phase becomes unreliable due to nesting effect. At the same time, there is a 
complex competition between SDW, CF phases, and <i-wave superconductiv- 
ity in two- and many-patch RG approaches. While the spin-flux (SF) phase 
is not stable at the mean-field level, it is identified as a possible ground state 
at J < in both RG approaches. With increasing t' the results of all three 
approaches merge: d-wave superconductivity at J > and ferromagnetism 
at J < become the leading instabilities. For large enough V the charge- 
density-wave (CDW) state occurs. 
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I. INTRODUCTION 



Despite many years of research the mysterious properties of underdoped high- 
temperature superconductors remain an unsolved and most challenging problem of strongly 
correlated electron physics. Recent theoretical proposals as well as experimental data trace 
the possible origin of the unconventional electronic properties of underdoped cuprates to a 
delicate competition between various types of ordering phenomena. Since metallic cuprates 
result from doping Mott-Hubbard antiferromagnetic insulators, it is natural that strong 
magnetic correlations persist for superconducting compositions but alternative ordering ten- 
dencies may be hidden as well. 

The relation between superconductivity and antiferromagnetism and the attempt to treat 
both types of order on the same footing was the basis of the phenomenological SO (5) theory 
as originally proposed by Zhang [1]. Yet, the observation of inhomogeneous charge and spin 
structures in static stripe patterns in rare-earth doped La 2 - x Sr :E Cu04 [2] and the appearance 
of an anisotropic rf-wave like pseudogap at temperatures well above the superconducting 
transition have pointed to the possibility of nearby instabilities towards other phases in 
competition or coexistence with rf-wave superconductivity. 

The relevance of stripe patterns was proposed early on (for a review see, e.g., Ref. [3]), 
and years later their existence was indeed experimentally verified in La 2 _ x _j / Ndj / Sr :E Cu04 [2]. 
These static stripe structures may in fact result from the lattice anisotropy connected to the 
Cu0 6 octahedral tilt pattern in the LTT phase of this compound [4]. Time-reversal symme- 
try breaking charge-flux (CF) [5-7] and spin-flux (SF, spin nematic) [8-10] phases were also 
proposed early on as possible ground states of the basic correlated electron Hamiltonians 
commonly used to model high-T c superconductors. 

The idea of circulating orbital currents (also called an orbital antiferromagnet state [5]) 
was recently revived as the underlying origin of the pseudogap phase in underdoped cuprates 
[11]. Two variants of orbital current structures are currently discussed: Varma [12] proposed 
a circulating current pattern in the Cu0 2 planes which preserves the discrete translational 
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symmetry of the lattice, while the more recently proposed <i-density wave (DDW) state 
[11,13] contains a staggered current pattern which doubles the unit cell. The latter state is 
in fact equivalent to the orbital antiferromagnet discussed a long time ago [6]. These time 
reversal symmetry breaking phases may furthermore coexist with rf-wave superconductivity. 
Yet, another instability connected to the spontaneous deformation of the Fermi surface 
(Pomeranchuk instability, PI) was recently proposed for the two-dimensional Hubbard model 
[14] . It appears tempting to suggest a possible relation between the Pomeranchuk instability 
and stripe pattern formation. 

Remarkably, the possible existence of orbital currents has most recently received exper- 
imental support. Angular-resolved photoemission with circularly polarized light identified 
intensity differences for left and right circularly polarized photons below the pseudogap 
temperature in Bi2212 [15] and c-axis oriented ordered magnetic moments were detected by 
spin-polarized neutron scattering in YBCO [16]. Both data sets can find a natural interpre- 
tation in terms of planar circulating current phases. 

The tendency towards the formation of charge flux phases is most likely beyond the 
physics of the standard Hubbard model. In order to allow for a larger variety of competing or 
coexisting phases it appears therefore demanding to explore some extensions of the Hubbard 
model. A natural choice is to consider the extended U-V-J Hubbard model (V is a nearest- 
neighbor interaction and J is the Heisenberg exchange coupling). This model was chosen 
early on to study orbital antiferromagnetism and spin nematic phases. In one dimension 
(ID), this model was studied within the bosonization technique [17,18]. However, in ID, 
one can not explore the possibility of charge- and spin-flux phases, but they were later on 
investigated in an extension to ladder systems [19]. In two dimensions, the U-V-J model 
was considered within mean-field approximation for J = [20,21] and for 1^ = [22]. The 
possible coexistence of antiferromagnetism and rf-wave superconductivity [23], and charge 
flux and rf-wave superconductivity [24] was discussed in this context. More recently, the 
possibility of coexisting phases was reanalyzed within a mean-field treatment of the extended 
Hubbard model which furthermore includes correlated hopping [25]. 
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Common to these studies is the need to go beyond the standard one-band Hubbard model 
to allow for a variety of alternative ordering phenomena. However, a mean-field treatment 
does not allow for an accurate treatment of fluctuations which may alter the complex ground- 
state phase diagram. Other analytical and numerical methods are needed for a reliable 
insight into the interplay of various types of ordering tendencies. One possible step forward 
in this direction has been the recent analysis of competing phases within dynamical mean- 
field theory [26]. 

Another possible route is to start from the weak-coupling regime, where magnetic or 
superconducting instabilities can result from nesting of the Fermi surface (FS) or van Hove 
singularities (vHS) in the density of states. The latter case is also relevant from a material 
oriented point of view since the FS of cuprates in an intermediate doping regime is close to 
vHS of the electron spectrum [27-29]. Schulz considered [30] different phases near van Hove 
band fillings within a mean-field-like analysis which, however, did not take into account the 
interplay of different electron scattering channels. The parquet approach was first applied to 
the extended Hubbard model with a nearest-neighbor hopping in Refs. [31-33], where besides 
the standard superconducting, spin- and charge- density- wave instabilities, the possibilities 
for spin- and charge-flux instabilities were also discussed. Recently, the standard Hubbard 
model (V — J — 0) was investigated within the renormalization-group (RG) approach on a 
patched FS using the Polchinsky flow equations [34], Wick-ordered RG equations [35], and 
a RG scheme for one-particle irreducible functions [36]. The possibility of the charge-flux 
phase as well as the Pomeranchuk instability for the standard Hubbard model was addressed 
in Ref . [37] . These latter methods have the advantage that they allow to take into account 
the contributions of the entire FS. At the same time, even the restriction of the momenta 
to the vicinity of van Hove points (the so called "two-patch" approach originally proposed 
in Refs. [38,39]) can capture already the essential physics and can give results which are in 
qualitative agreement with the RG approaches on patched Fermi surfaces, see Ref. [40]. This 
two-patch approach allows to include the contribution of particle-hole scattering at small 
momenta which is expected to be important for the determination of the complete ground- 
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state phase diagram. The role of this type of scattering was also explored within the RG 
approach with a temperature cutoff on the patched FS [41]. Although the phase diagram of 
the U-V-J Hubbard model was investigated previously within the two-patch RG approach, 
only the cases of J = [21] and arbitrary J with nearest-neighbor hopping dispersion [42], 
but without the contribution of the small-momenta particle-hole scattering were considered. 
The same model within the many-patch approach was not systematically studied so far. 

In this paper we reconsider the phase diagram of the extended U-V-J Hubbard model 
with special emphasis on the possibility of flux phases. We employ the mean-field approxi- 
mation, two- and many-patch RG approaches. The paper is organized as follows: In Sect. 
II we start with the mean-field phase diagram of the U-V-J Hubbard model. In Sect. Ill we 
describe the two- and many-patch RG schemes and compare the results obtained with both 
approaches to the mean-field phase diagram. In Sect. IV we complement the RG analysis 
by symmetry arguments. Finally, in Sect. V we summarize the main results and conclude. 

II. MEAN-FIELD ANALYSIS 

We consider the extended U-V-J Hubbard model on a square lattice as given by 

ija i (ij) (ij) 

where the hopping amplitude ty = t for nearest neighbor sites i and j and = — £' for 
next-nearest neighbor sites (£,£' > 0); N e is the total number of electrons. U > is the 
on-site and V the nearest-neighbor interaction, J is the Heisenberg exchange coupling, and 

n i = ^ j c lg C i<Ji Sj = — ^ ] Cj a (^aa'Cia'] 

(T ua i are the Pauli matrices. We introduce operators which correspond to different types of 
order 

6 C DW« = ^(-l)^(cU CT -n/2), 
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Of« = ^<t( C U-1/2), (2) 

where A^- = +1 (—1) for j — i + 8 x (5 y ) and zero otherwise; denote unit vectors which 
connect to nearest neighbor sites in x and y directions respectively. (—1)* = (—l) tx+lv , z = 4 
is the number of nearest neighbor sites on the square lattice, and n is the electron density. 
Nonzero average values of the operators Ocdw(^) an d Osdw(*) correspond to charge- and 
spin density wave states, Odsc(0 to rf-wave superconductivity, Oqf and Osf to states with 
charge flux (the orbital antiferromagnet state) [5,6,13] or spin flux (spin nematic state) 
[8-10], respectively, and finally Of corresponds to a ferromagnetic state. We decouple the 
interactions by introducing the order parameters O m = (O m ) (m = CDW, SDW,... ), and 
arrive at the mean-field Hamiltonian 

H = £ kcLc k(7 + A m O m (k) + ^ (3) 

where 

e k = — 2t(cos k x + cos k y ) + At' cos k x cos k y — /i. (4) 
O m (k) are the Fourier transforms of the operators in Eq. (2), A m = O rn T m and 
r CDW = 8V-U; r SDW = U + 2 J; r dSC = 3 J - W; 

(5) 

r CF = 3 j + 4V] r SF = 4\/ - j ; r F = u - 2 j. 

At finite temperature T the self-consistency conditions lead to the mean-field equations 
(see e.g. Ref. [8]) 
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2N E* 2T ' 

k,a=± ka 

n = 1 V tanh^, (6) 

k,a=± 

where E' denotes the summation over momenta in the magnetic Brillouin zone defined by 
cosk x + cos k y > (the lattice constant has been set to unity); </> k = 1 for SDW and CDW 
phases and k = / k = (cosk x — cosk y )/2 for CF, SF, and dSC phases, 

ez = £ -±±^±* + f^ k - e k+Q r + m^ m m 

for m = SDW, CDW, CF or SF and 

EE? = \ V (^-^ +Q ) 2 + 4^ sc . (8) 

Q = (ir, 7r) is the wavevector for staggered order. For the ferromagnetic phase the mean-field 
equations read 



ka 



k,«=± ka 

El a = £k - «A F (9) 



where / ka = f(E^ a ) is the Fermi function. Non-trivial mean- field solutions exist only for 
T m > 0, which therefore determines the borders of absolute stability of the corresponding 
phases in the phase diagram. 

We also allow for coexistence of antiferromagnetism, ci-wave superconductivity, and 
charge-flux phases. The corresponding MF equations are the generalization of those for 
the coexistence of SDW and dSC phases [23] and CF and dSC phases [24], and have the 
form 

r SDW 



2N 

k,a=± 



Y -^tanh%fl + £k + £k+Q V 
^ £ ka 2T V 9ka J 



1 = ^T y #tanh%, 
2N E ka 2T 

k.,a=± 



n = 



1 




(10) 



where the electronic spectrum is given by 



E** = yj(ei + 4 +Q )/2 + ^sdw + tii^CF + 



dSC 



) + (£k + £k+q)^fca/2, 



^ka = ay/ (e k - e k+Q ) 2 + 4(A 2 DW + f£A 2 dSC ). 



(11) 



We note that a nonzero expectation value of both, antiferromagnetic and superconducting 
order parameters, produces also a nonzero 7r-triplet pairing amplitude with momentum Q 
(see e.g. Ref. [23]). Analogously nonzero values of charge flux and d-wave superconducting 
order parameters lead to a finite singlet pairing amplitude with momentum Q (so called 
77-pairing). Here, we do not include ir- and 77-pairing in the mean-field analysis, since it was 
verified numerically [23] that the corresponding amplitudes are significantly smaller than 
the amplitudes of the corresponding "parent" order parameters. 

The mean-field equations (6), (9), and (10) were solved numerically for T — 0, U = 4t, 
t' = O.li, and band filling n = 0.92 (which is the van Hove filling for this value of t'/t, 
see the next section). Since the regions of stability of different phases overlap, we compare 
the energies E = (H) + jiN e to identify the mean-field ground state. The resulting phase 
diagram is shown in Fig. 1. At positive J and not too large V we find coexistence of a spin 
density wave with d-wave superconductivity. However, the fraction of superconductivity in 
this phase is quite small (for V = 0, J = 0.5U we have Osdw = 0.40 and Odsc = 0.06). 
With decreasing V the superconducting fraction increases and becomes substantial close to 
the transition into the pure dSC state (for V = — 3C//4, J = — 0.05U we have Osdw = 0.18 
and Odsc — 0.12), and finally at the boundary between the dSC and SDW+dSC phases the 
antiferromagnetic order parameter vanishes continuously. Therefore, this transition as well 
as the transition from the SDW into the SDW+dSC phase is of second order. All other phase 
boundary lines in Fig. 1 correspond to first order transitions. The CF, SF, and CF+rfSC 
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phases are not energetically favorable at all on the mean-field level, the dSC phase is stable 
only for V < —2t. The mean- field phase diagram at other fillings or other t'/t is similar. 
In the special case t' — and half-filling, the coexistence region of the SDW+dSC phase is 
replaced by a pure SDW state. 

Away from half-filling the phases with ordering wavevector Q = (ir, ir), i.e. CDW, SDW, 
CF, and SF phases are usually unstable towards phase separation (see e.g. the discussion 
in Refs. [43,44]). To analyze this possibility, we calculate the isothermal compressibility k t 
defined by 



A negative k t necessarily implies the thermodynamic instability of the considered homo- 
geneous phase and a possible tendency towards phase separation. The numerical analysis 
shows that all above phases are indeed unstable except for the SDW+dSC phase in the 
region J <C U and V — > — U where the fraction of the SDW order parameter is small. 

The above mean-field results provide a rough qualitative picture for the regions of dif- 
ferent stable phases in the ground-state phase diagram of the U-V-J Hubbard model. More 
detailed information is extracted from the weak-coupling RG treatments discussed below. 



The tight-binding spectrum (4) leads to vHS in the density of states arising from the 
contributions of the points = (vr, 0) and k B = (0,7r). These singularities lie at the FS 
if n — — At'. For t' — the corresponding filling is uvh = 1 and the FS is nested, but the 
nesting is removed for t'/t > 0, when nyu < 1- The shape of the FS at different t'/t and 
van Hove band fillings is shown in Fig. 2. 

The two-patch approach [38-40] considers the fillings which are close to the van Hove 
band fillings. At these fillings the density of states at the Fermi energy and the electron- 
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III. RENORMALIZATION-GROUP ANALYSIS 



A. Two-patch analysis 
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electron interaction vertices at momenta k = \s.a,b contain logarithmical divergencies coming 
from the momentum integrations in the vicinity of the points k = \s.a,Bi an d therefore one 
can expect that the contributions from the vicinities of these points are the most important 
for the calculation of the renormalized electron-electron interaction vertices. 
We start by rewriting the U-V-J Hamiltonian in momentum space 

H = Y £ kcLc k a + Y Y £( kl > k2 ' k3 ' k 4) C I 1 a C k 2CT ' C k3^ C k4 CT ^k 1 +k 2 -k 3 -k 4 (13) 

where 



2N 2 

k<7 kik2k;ik4 (Tit' 



Ckcr = 



g(k u k 2 , k 3 , ki) = U + (V- J/4) 7 k 2 -k 3 - (^/2)7k 3 - kl (14) 

7k = 2(cos/c a: + cos k y ) and the Kronecker ^-symbol ensures momentum conservation. Since 
we restrict the momenta to the vicinity of k^ and kg, it is convenient to introduce new 
electron operators and 6k by 

a k ~ kA ,a keO(A) 
6 k -k BlCT k e 0(B) 

here 0(A) = {k : |k — k^| < A} and similar for B; A is a momentum cutoff parameter. We 
expand the spectrum near the van Hove points 

£k4+ P = £p = -2t(sin 2 ippl - cos 2 ipp 2 y )-Jl (15a) 
= -2tp + p_ - Jl 

£k B+P = el = 2t(cos 2 pp 2 x - sin 2 pp 2 y )-Ji (15b) 
= 2tp + p^ - Jl 

where cos(2<^) = R = 2t'/t, /! = // + At', p± = p x sin ip ± p y cos if, p± = p x sin ip ± p y cos tp 
and rewrite the Hamiltonian in the form 

h = Y £ i a U a p° + Y e ? b U b p° 

pa pa 

2-kH 

+ _ /V2~ Y \-9 l ^) a pia h p 2 a< a V3a' h Pi° + 92 (A)4>l^p 2 <7' & P3^'<W](Wp2-P3-P4 
Pi, aa' 

+ ^4(A)a Pia a P2CT ,ap 3 a'ap 4CT + a 6]<J Pl 

+P2-P3— P4 

(16) 
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where e p ' B = £k AB + P , A = In (A/ max(p i+ , pi_,p i+ ,pi_, \fJ>\/t); the summation is restricted to 
momenta |pj| < A. Neglecting the weak (non-logarithmical) dependence of the vertices gi 
on momenta at |pj| < A, we obtain the bare values of the vertices g® which are independent 
of A and given by 

gl = g(k A , k B , k A , k B ) = g (l - 4V/U - J/U), 
g% = g(k A , k B , k B , k A ) = g (l + 4V/U + J/U), 
g° 3 = g(k A , k A , k B , k B ) = <?o(l - W/U + 3J/U), 

g° 4 = g(k A , k A , k A , k A ) = g (l + AV/U - 3J/U). (17) 

where g = U/ (AttH) is the dimensionless coupling parameter. Note that at half filling and 
t' = the contributions of the flat, nested parts of the FS (which should be considered 
separately, see, e.g. Ref. [33] and references therein) is neglected. The contribution of these 
nested parts of the Fermi surface can be considered within the many-patch approach, as we 
discuss in the next section. 

To obtain the dependence of the vertices gi on A we apply the RG analysis. We start with 
the bare values of the vertices at momenta far from van Hove singularities, i.e. with g®, and 
integrate out the fermions a p with momenta Ae _A < p± < Ae~ A_dA and fermions b p with 
momenta Ae~ A < p± < Ae~ x ~ dX at each RG step. We consider first the one-loop corrections 
which contain particle-hole (ph) and particle-particle (pp) bubbles at small momenta and 
momenta close to Q. The results for these bubbles can be summarized as follows: 

™ y 1 ~ ffffl ~/( £ p+ q ) = A 2 (18a) 
< +Q = E /(£ f]-y^ ^^min ( A 2 ,2, Q ), (18b) 

yrpp _ 1 ~ f( E p) ~ /( £ p+q) CQ 

q+Q " p — — "2^' m 

Ae- A <p±<A p p+q 



where 



c = l/sin(2p) = 1/Vl -R 2 , 
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c Q = tan _1 (i2/Vl - R 2 )/R, 

zq = ln[(l + VT^TP)/R}. (19) 

The contribution of the slices Ae~ A < p±,p± < Ae~ A_dA is obtained by taking the derivatives 
dUi/dX. Since the contribution to the particle-hole (Jl^ 1 ) channel is concentrated near the 
FS, we use the weaker cutoff condition Ae~ A < p + < A or Ae~ A < p_ < A for this channel 
which gives 

yiph _ / ( £ P ) ~ / ( £ P+q) _ Z 0^ (on \ 

q " x ^ > ^P " £ P+Q _ 2 ^ 

Ae~ A <p + <A or Ae~ A <p_<A p 

with 2 = c . A more accurate treatment of the ph channel requires the parquet summation of 
the one-loop diagrams, it was shown, however, in Ref. [40] that the corresponding corrections 
do not change qualitatively the results. Equation (20) can be justified also by considering 
the temperature flow of the vertices on a patched FS as proposed in Ref. [41] (see also the 
next section). 

While the bubbles Hq (at A < 2zq) and n"o P contain squared logarithms, the bub- 
bles Uq (at A > 2zq), I1q /i and I1q contain only single-logarithmical divergencies. Strictly 
speaking, the RG approach does not perform a correct summation of these subleading diver- 
gencies. However, one can not simply neglect the single-logarithmical contributions (n^q 
at A < 2z Q , Uf and nf ) in comparison to the squared-logarithmical ones, since it amounts 
to ignore some electron scattering channels and therefore would exclude the possibility of 
related orders, see Ref. [40] for a discussion. Mathematically, this difficulty manifests itself 
in the corresponding phases by the growing of the vertices which are multiplied by single 
logarithms, whose contribution can therefore become comparable to the squared- logarithmic 
contributions [40]. 

As we will see below from the solution of the RG equations, in the weak-coupling regime 
at small t'/t 1 the contribution of the subleading corrections Hq 1 and I1q is small so 
that they do not change substantially the flow of the coupling constants. With increasing 
coupling, the contribution of single-logarithmical terms is comparable to the contribution 
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of the leading squared-logarithmical terms which reflects the possibility of a corresponding 
ordering phenomena (e.g. ferromagnetism) in the strong-coupling regime, although the latter 
can not be explored within the weak-coupling RG approach. 

At intermediate t'/t the situation changes. In this case the contributions of different 
channels (containing single- and double-logarithmical terms) are comparable already in the 
weak-coupling regime. We have checked however, that even for intermediate t'/t the quali- 
tative results of the two-patch RG approach do not depend on the cutoff procedure, which 
means that the treatment of the subleading single-logarithmical terms within the same RG 
procedure for intermediate t'/t is qualitatively reliable. Therefore, we take into account the 
contribution of both, single- and squared-logarithmical terms in Eqs. (18) and (20). At the 
same time, the contribution of the subleading single-logarithmical term in Eq. (18a) can be 
safely neglected, since it is multiplied by the same vertex as the leading one and therefore 
is always subleading. 

We determine the RG equations for the vertices gi(X) in the form [38-40] 

dgi/dX = 2d 1 (X)g 1 (g 2 - gi) + 2d 2 gig4 -2d 3 g 1 g 2 , 
dg 2 /dX = di(A)(#2 + ol) + 24(5-1 - #2)54 - d 3 (gf + g%), 
dg 3 /dX = -2d (X)g 3 g 4 + 2d 1 (X)g 3 (2g 2 - g x ), 

dg A /dX = -d (X)(gl + g\) + d 2 {g\ + 2 9l g 2 - 2g\ + g\), (21) 

where 

d (X) = 2c A, d 2 = 2z ; d 3 = 2c Q ; 

di(A) = 2min(A,z Q ). (22) 

Eqs. (21) have to be solved with the initial conditions <?i(l) = §f . Note that we neglect here 
the renormalization of coupling constants which arise from non-logarithmical contributions 
at A < 1, this renormalization is small provided that the condition g° <C 1 is satisfied. 

Eqs. (21) coincide with the temperature cutoff RG approach for one-particle irreducible 
functions on the patched FS [41] (see Sect. IIIB) in the approximation that the vertex is 
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constant (i.e. does not depend on the patch) for the momenta in the vicinity of the vHS 
(for | pj | < A) and zero far from it (when |pj| > A at least for one i = 1...4). Therefore, the 
real dependence of the vertices on the momenta along the FS is replaced by a step function 
in the two-patch approach. At the same time, the momentum dependence of the vertices 
in the direction perpendicular to the FS is treated correctly through the scaling variable A. 
Note however, that the momentum dependence of the electronic spectrum within each patch 
is correctly taken into account in the two-patch approach. 

Strictly speaking, some of the two-loop diagrams (e.g. the self-energy corrections) give 
contributions of the same order as the single-logarithmical terms which we treat here and 
they should be taken into account as well. However, the analysis of such contributions is 
rather involved and the subject for future work [45]. Note that this difficulty is also "hidden" 
in the RG approaches on a patched FS where additional logarithmical divergencies arise from 
the momentum integration along the FS. 

In order to explore the possible instabilities of the system, we consider the behavior of 
the zero-frequency, time-ordered response functions at zero temperature 



where O m (q, r) denote the Fourier components of the operators (2) in the Heisenberg rep- 
resentation. Besides the operators in Eq. (2) we also test, following Refs. [21,42], the 
susceptibilities which correspond to the operators 



oc 




(23) 



— oo 





(24) 
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where are the nearest neighbor sites of site i. The susceptibility Xq(0) is related to the 
isothermal compressibility k t by 

«r = XQ(0)/n 2 . (25) 

The bond-charge order parameter r characterizes an instability towards a spontaneous de- 
formation of the FS [42], i.e. anisotropy in x and y directions, the so called Pomeranchuk 
instability [14]. Analogously, a bond-spin order parameter A characterizes a spin-dependent 
Pomeranchuk instability. The operators n and rj correspond to triplet and singlet pairing 
with momentum Q = (7r,7r). 

Picking up the logarithmical divergences in A we obtain the RG equations for the sus- 
ceptibilities in the same approximations as discussed above (cf. Refs. [30,38-40]) 

d Xm (A)/dA = 4 m (A)T^(A), (26) 
dlnT m (A)/dA = 4 m (A)r m (A), 

where the coefficients Y m are given by 

Tcdw = 92—93 — 2#i; Tsdw = 92 + 93', 
r C F = 92 + 93 - 2#i ; rsF = fi'2-fi'3; 
I\ = 9i - 92] = -g 1 - g 2 ] 

(27) 

r Q = 9\ - 2# 2 - 9a] T F = g 1 + g 4 ; 
T T = -91 + 2g 2 - g4] F A = g 4 - gi, 
Tdsc =93—94- 

In Eqs. (26) a m = for the dSC phase; a m = 1 for SDW, CDW, SF, and CF phases, a m = 2 
for F, Q, r and A phases, and a m = 3 for n and rj phases. Eqs. (26) have to be solved with 
the initial conditions T m (0) = 1, Xm(0) = 0. At A = the vertices (27) coincide with those 
considered in the mean-field approach in section II (see Eq. (5)). 

Numerical solutions of Eqs. (21) show, that at a critical value A c of the scaling parameter 
A some of the vertices and susceptibilities are divergent. We analyze the behavior of the 
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coupling constants g\ to g± for A — > A c representing it in the form (fa. ..fa) where fa — 0, + 
or — describes the behavior of the coupling constant gf the "plus" ("minus") sign means 
relevant in the RG sense and tends to +00 (—00) and zero means irrelevant. To identify the 
leading instabilities we calculate the inverse susceptibilities for A < A* (A* is the value of the 
scaling parameter A when the largest absolute value of the coupling constants \gi\ exceeds 
unity) and then use a linear extrapolation of the inverse susceptibility (see Fig. 7 below). 
From the RG point of view, at A ~ A* the "one-dimensional-like" behavior of the coupling 
constants considered above changes to two-dimensional and also the parts of the FS far from 
the van Hove points become non-negligible. However, since the values of the susceptibilities 
are large already at A*, the corresponding critical region is narrow. 

For a given A c the size A of the patches is restricted by A <C n and ln(4/A) <C A c . 
The latter criterion follows from the condition that the contribution of the electrons with 
\k±\ < A to particle-hole and particle-particle bubbles is dominant (see, e.g. Ref. [42]). We 
choose A = 1 and require A c 3> In 4 ~ 1. Since A c decreases with increasing interaction, this 
criterium restricts the values of the interactions where the two-patch RG approach is valid. 

B. Many-patch analysis 

In the many-patch analysis we follow the temperature-cutoff renormalization group 
for one-particle irreducible (1PI) Green functions proposed recently by Honerkamp and 
Salmhofer in Ref. [41]. This version of the RG uses the temperature as a natural cutoff 
parameter, allowing to account for both the excitations with momenta close to the Fermi 
surface and far from it, which is necessary for the description of the particle-hole instabilities 
with zero momentum transfer, e.g. ferromagnetism, phase separation, and the Pomeranchuk 
instability. Neglecting the frequency dependence of the vertices, which is justified in the 
weak-coupling regime, the RG differential equation for the temperature- and momentum- 
dependent electron-electron interaction vertex (see the diagrammatic representation in Fig. 
3), has the form 
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-^V T (ki, k 2 , k 3 ) = - VrQzi, k 2 , p)L pp (p, -p + ki + k 2 )F T (p, -p + k x + k 2 , k 3 ) 

p 

+ [ - 2^r(ki, p, k 3 )y T (p + ^ - k 3 , k 2 , p) + V^ r (ki, p, k 3 )F T (k 2 , p + k x - k 3 , p) 

k 

+Vr(ki, p, p + ki - k 3 )V r T (p + ki - k 3 , k 2 , p)]L ph (p, p + ki - k 3 ) 

+ ^(ki, p + k 2 - kg, p)L ph (p, p + k 2 - k 3 )Vr(p, k 2 , k 3 ) (28) 

k 

where 



V(k,k') 



£k — £k' 



L pp (k,kO = ^ (£k) + ^ (£k,) (29) 
£k + £k' 

and /t( £ ) = df(e)/dT. Eq. (28) has to be solved with the initial condition Vr (ki, k 2 , k 3 ) = 
g(k 1 , k 2 , k 3 , k x + k 2 — k 3 ) where the initial temperature T is of the order of the bandwidth. 
The evolution of the vertices with decreasing temperature determines the temperature de- 
pendence of the susceptibilities according to [41,37] 
d 



XmT = ^T mT (k')T mT ( T k' + Q m )L pPiP/l (k', T k' + q m ), (30) 
k' 

-^T mT (k) = =F ^ T mT (k')r mT (k, k')L pp , pA (k', T k' + q m ) 



dT 

df 



where 



V T (k, k', k' + q m ) - 2V T (k, k', k + q m ) for m = CDW, CF, Q and r, 
()l r(k, k') = { \/ T (k, k', k' + q m ) for m = SDW,SF,F, and A, (31) 

V T (k, — k + q m , k') for m = 7r, 77, and dSC. 

q m = Q for CDW, SDW, CF, SF, ir and 77, and q m = otherwise. The upper signs in Eq. 
(30) refer to the particle-particle response (ir, 77 and dSC), the lower signs to the particle- hole 
response. The initial conditions for Eqs. (30) are 

{cosA^ — cosk v for CF, SF, PI, A, dSC, and n, 
(32) 
1 otherwise, 

and Xm,T = 0. To solve numerically Eqs. (28) and (30), we use the discretization of mo- 
mentum space in 32 patches and the same patching scheme as proposed in Ref. [41] (we 
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have checked in selected cases that increasing the number of patches to 48 does not change 
the results). With account of the symmetries of the square lattice, this reduces the above 
integro-differential equations to a set of 1920 differential equations which were solved nu- 
merically. We use the value of the starting temperature T = 12t, which is slightly larger 
than the bandwidth, and stop the flow of the coupling constants where the largest cou- 
pling constant V max = 18t. As for the two-patch analysis, we extrapolate the resulting 
inverse susceptibilities and calculate the critical values of the scaling parameter A™ where 
the extrapolated inverse susceptibilities vanish. Due to this procedure, the results of RG 
analysis do not depend strongly on V max . Note that the initial /c-dependence of the response 
functions (32) is slightly changed during the renormalization-group flow: responses with 
d-wave symmetry (CF, SF, PI, A, dSC, and n) acquire g-wave and higher-order harmonics, 
while responses with s-wave symmetry acquire additional extended s-wave (cos k x + cos k y ) . 
However, these additional corrections are small. 



C. Results of the RG analysis 

For the calculations we choose the interaction strength U — 2t since for stronger inter- 
actions and moderate \J\/U and |V|/C/ the RG approach becomes unreliable. The phase 
diagrams calculated with this value of U, different values of t'/t = 0, 0.1, 0.3 and the 
corresponding van Hove fillings are presented in Figs. 4-6. The phase boundaries, deter- 
mined from the two-patch approach are shown by thin solid and dashed lines. Solid lines 
separate the phases with different behavior of the coupling constants, while dashed lines 
separate phases with the same behavior of the coupling constants and are determined from 
the condition of the equality of the corresponding critical values A = A™ where the linearly 
extrapolated susceptibilities vanish. 

We consider first the results of the two-patch approach. Since this approach is applicable 
for not too small A c (see Section IIIA), we consider only the region of the phase diagrams 
with A c > 2, which is bounded by bold lines. In the most part of the phase diagrams the 



18 



typical values of A c range from 3 to 5. The dependencies of the inverse susceptibilities on 
the scaling parameter for selected parameter values are shown in Fig. 7. For not too large 
values of \J\ (J < 0) and for t' — or t' — O.lt (Figs. 4,5) the two-patch approach predicts 

that the spin-flux phase of the type (0 H ) is the leading instability. With increasing 

V this phase is replaced by a phase which has comparable CDW and CF susceptibilities 

(CDW phase) and the behavior of the coupling constants ( 0+). With increasing J up 

to J > 0.25£7 the CDW susceptibility dominates but also Xr and xsf are large (we term 

this phase CDW'). The corresponding flow of the coupling constants is (OH ). At large 

negative V a phase of type (0 — 0—) occurs where the largest diverging susceptibility is xq 
which, according to Eq. (25), implies a divergent isothermal compressibility k t . A divergence 
of kt can be attributed generally to different phenomena: a metal-insulator transition [46] 
or phase separation. In the considered parameter region the most natural explanation is 
phase separation induced by a large negative V. For large enough \ J\ (J < 0) we obtain 
the ferromagnetic state (F) with the coupling constant flow (+ + 0+); for large positive V 
Xa is the most diverging susceptibility. However, for small t' < 0.3t these instabilities are 
outside of the weak-coupling region of the phase diagram. The dashed area corresponds 
to the frustrated regime where the critical value of the scaling parameter A c > 20 and the 
spin-density-wave instability is strongly suppressed by negative J. Such behavior arises from 
the competition between antiferromagnetic and superconducting fluctuations on one side 
and ferromagnetic fluctuations on the other side. As a result different types of instabilities 
almost "cancel" each other. 

Four other phases (SDW, dSC, CF and PI) have the same flow of the coupling constants 
(0 + H — ). The magnitudes of the susceptibilities Xdsc, Xsdw, Xcf in the region J > are 
close to each other (see Fig. 7b), which naturally implies a close competition between these 
states. With increasing J we find a PI phase where the susceptibility Xr is largest. However, 
this instability also always appears outside of the weak-coupling region of the phase diagram 
where the equations (21) are valid. 

As in the mean-field solution, away from half filling one expects phase separation or 
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the formation of inhomogeneous structures of all ordered phases with wavevector Q (CDW, 
SDW, CF and SF). Note however, that this filling- induced type of phase separation should 
be contrasted with the interaction-induced phase separation in the PS region of the phase 
diagram, which is present even at half-filling (t' = 0) and may not be magnetically ordered 
or superconducting. 

The results of the many-patch approach are shown by different symbols, explained in 
the figure captions. Solid symbols correspond to different types of ordering tendencies. 
Open symbols correspond to frustrated behavior, where the many-patch RG could not reach 
Knax = 18t for At = ln(£/T)/2 < 4. Boundaries between some phases (e.g. CDW and SF, 
dSC and PS for J > 0) almost coincide in the two- and many-patch approaches. At the 
same time, the many-patch approach predicts a much broader region of stability for the SDW 
phase and an almost vanishing region of stability for the CF phase. Therefore, for t' — 
and t' — 0.1 taking into account more patches on the FS (and therefore the contribution of 
nested parts) moves the predicted phase boundaries closer to those of the mean-field results. 
At the same time, the charge flux, as well as the <i-wave superconducting response in the 
many-patch approach are substantial within the antiferromagnetic phase, which implies a 
close competition between the above ordering tendencies. For a rough estimate, where the 
charge response is closest to the antiferromagnetic one, we bound the region with xcf > 
(2/3)xsdw by thin crosses. 

At J < many-patch approach does confirm the possibility of a spin flux phase in almost 
the same parameter region as determined from the two-patch approach. With decreasing 
V this spin flux phase is however replaced by a (partially frustrated) ferromagnetic phase, 
which fills most of the region where the frustrated SDW order was expected in the two-patch 
approach. The presence of the ferromagnetic phase at J < resembles closeby the mean- 
field prediction in this parameter range. The d-wave superconducting region for J < is 
replaced by phase separation in the many-patch approach. 

Note that with decreasing of U the region of <i-wave superconductivity at J > grows, 
while increasing of U favors the SDW instability at J > 0. Similarly, at J < decreasing U 
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shifts the balance between SF phase and ferromagnetism towards the spin-flux phase. Other 
phases (e.g. CDW and PS) are essentially not influenced by varying U. 

When increasing t' up to t' = 0.3t the phase diagram changes (see Fig. 6). In the two- 
patch approach, at J > 0.05U d-wave superconductivity becomes the leading instability. 
For large positive V we obtain either CDW and A instabilities (CDW phase on Fig. 6) or 
CDW and PI instabilities (CDW phase) depending on J. For large negative V we again 
find interaction-induced phase separation. Most of the region with J < 0.05U and moderate 
|V| is frustrated (A c > 20), which is again the result of a strong competition of different 
ordering tendencies. 

The many-patch approach also gives the strongest tendency to rf-wave superconductivity 
at J > 0.05C/, but the tendency towards ferromagnetic order in a large part of the frustrated 
region J < 0.05U. Similar to t' — and t' = O.lt, at t' = O.St, the rf-wave superconducting 
region for J < predicted by the two-patch approach is replaced by phase separation in 
the many-patch approach. The boundary of the CDW phase also almost coincides for the 
many-patch and two-patch approach. Therefore, at t' = 0.3t the predictions of the two-patch 
approach are closer to the results of the many-patch approach, than for smaller t' . This is 
connected with the absence of nested parts of the Fermi surface, as discussed above. 

In summary, the RG analysis gives a much richer phase diagram than anticipated from 
the mean-field results of the previous section. At the same time, for small t' the results of the 
many-patch RG approach, which takes into account the contributions of the whole FS and 
therefore treats the renormalization of the interactions in a more accurate way, are closer 
to the mean-field predictions than the results of the two-patch approach. At intermediate t' 
the two-patch approach becomes more reliable and the results of both approaches are close. 
In the next section we supplement the above RG studies by an 5*0(8) symmetry analysis. 
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IV. SYMMETRIES AND ASYMPTOTIC BEHAVIOR 



Additional insight to the phase diagram is obtained by considering the symmetries of 
the Hamiltonian (1) with respect to the 5*0(8) symmetry group of transformations of the 
4-fermion states |k f ), |k j), |k + Q f) and |k + Q j) (see e.g. Ref. [49]). This group is 
generated by the so(8) algebra of 28 operators which are given in part by 

o m = ^J2 d ^ ( 33 ) 

i 

with O m (i) as defined above in Eqs. (2) and (24). For a full list of operators of the so(8) 
algebra we refer to Ref. [49]. From the operators O m of the so(8) algebra we construct the 
rotation operators 

Rm(a) = exp(iaO^ + iaO m ) (34) 

where a is a real number. From the symmetry point of view, the operators O m = n, rj, A, 
and f are most useful. As discussed in Refs. [49,21], the operator R n performs a rotation 
between SDW and dSC states (which is the basis of the SO (5) theory [1]), the operator R v 
between dSC and CF states, the operator Ra between SDW and CF phases as well as CDW 
and SF phases while the operator R T rotates between SDW and SF phases and between 
CDW and CF phases, respectively. 

The symmetry of the Hamiltonian (1) under these operations is as follows. The non- 
interacting part of the Hamiltonian is invariant (H = R m H R^) under the rotations Ra 
and R T and in the nesting case = — £k+Q under R n and R v . The interaction is invariant 
only under the rotations R v for V = J = 0, i.e. the Hubbard model, which was originally 
discovered by Yang and Zhang [50] but the interaction is not invariant with respect to the 
other rotation operations. 

The symmetry of the restricted two-patch Hamiltonian (16) is considerably higher [42]. 
For this restricted Hamiltonian the interaction part, and therefore the entire Hamiltonian 
at half-filling and zero t', is invariant with respect to the operations R w , R v , Ra or R T on 
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special lines in the V-J plane (at fixed U). The corresponding lines are supplied by captions 
in Fig. 4, corresponding to the type of the symmetry. The symmetry lines exactly coincide 
with the boundaries between different phases in the two-patch approach, since the corre- 
sponding susceptibilities Xm(A) in this approach are identical on these lines for arbitrary A. 
At the same time, as we have seen in Section III, most of these symmetries (those which are 
connected with the SDW order) are broken by the contribution of the nested parts of FS. 
Away from half-filling or at finite t' the non-interacting part of the Hamiltonian becomes 
invariant only with respect to the rotations Ra and R T . 

To get more insight into the posibility of different ordering tendencies, we rewrite the 
Hamiltonian in terms of the operators of the so(8) algebra. We restrict ourselves to the 
contribution of the vH points, since the corresponding analysis for the whole FS becomes 
very complicated and requires numerical diagonalization of the corresponding Hamiltonian. 
Taking into account that the non-interacting part of the Hamiltonian (16) for vH points, i.e. 

vanishes at the vH filling (where e^ A = = 0), the Hamiltonian is rewritten as 

Heff = 9i(0 2 CBW - 2 CF ) + g 2 (Q 2 - r 2 ) + g 3 (d 2 CDW + 2 CF ) (36) 
+gi(Q 2 + T 2 )-2Q(g 2 + g A ) 

where Q = Q+l/2 and g 1 ... g 4 are the vertices determined in the two-patch approach (Sect. 
IIIA). These vertices can be also considered as those obtained from the many-patch analysis 
for the corresponding momenta, provided that the contribution of electrons with momenta 
far from the vH singularities is small, i.e. t'/t is not close to zero. Under these conditions, 
the Hamiltonian (36) can be considered as an effective Hamiltonian of the RG procedure in 
the previous section, which acts on the space of 16 states \sa, sb) where sa,b = 0, T, I or f j 
denotes the electron states with vH momenta or k B . 

First we consider the eigenstates \m) of the operators O m , which are easily expressed in 
terms of the states \sa, sb) and we obtain the expectation values of the Hamiltonian in the 
corresponding states as 
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E m {\) 



{m\H eff \m) 
(m\m) 



C(A)-r m (A)/2 



(37) 



where C(A) = — (gi + 2g 2 — 2g 3 + g^/2 is independent of m and T m (A) is defined in Eq. 
(27). Therefore the states which were identified in the two-patch approach as having the 
largest susceptibility (largest T m (A) for A — > A c ) have also the lowest effective energy among 
the eigenstates \m). 

For m ^ F or Q the eigenstates \m) of the operators O m are not the eigenstates of 
the effective Hamiltonian (36), since the operators O m generally do not commute with this 
Hamiltonian. To find the eigenstates of the Hamiltonian, we represent it as a matrix with 
respect to the states \sa,sb) and diagonalize it. In fact, the subspaces with even (odd) 
number of particles 



are not mixed by the operators of the so(8) algebra and we can consider them separately. In 
the following, we focus on the even subspace where the interaction part of the Hamiltonian is 
non-trivial. The resulting energy levels and the eigenfunctions are presented in Table 1. The 
states \Vi) are the triplet of eigenvectors of the operator Of with eigenvalues 1,0 and — 1, 
respectively. From the states \Ui) one can form a doublet of vectors \Ui) ± \U 2 ) which are the 
eigenvectors of the operator Q with eigenvalues =Fl. The other states are not eigenvectors of 
any operator in Eqs. (2) and (24). Considering the flow of the energy levels of the respective 
states which are given as a function of the coupling constants g t in the second column of 
Table 1, we obtain the correspondence between the lowest energy eigenstates of the effective 
Hamiltonian in the two-patch RG procedure and the phases identified in Sect. III. We find 
that in the F and SF phases the triplet \Vj) has the lowest energy. It is natural to associate 
the states (Vi^), which have eigenvalues ±1 of the operator Of, with the ferromagnetic phase 
and the state \ V 2 ) with the SF phase. Similarly, in the PS phase the states |C/j), in the CDW 
and A phases the state |Ci), and in the CDW phase the state |C 2 ) have the lowest energies. 
Finally, the \W) state is common for SDW, CF, dSC, and PI instabilities, i.e. phases which 



(38) 
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have the coupling constants flow of (0-H — ) type. The resulting correspondence is presented 
in the right column of Table 1. 

To clarify possible types of orders which correspond to the eigenstates \k) of the effec- 
tive Hamiltonian (36), we expand these eigenstates in terms of the eigenvectors \m) of the 
operators O m by calculating the scalar products C mk = (k\m) listed in Table 2. Generally 
the states \k) are a mixture of the states which correspond to different order parameters. 
So, the state \C±) in Table 1, which was identified above with CDW and A phases, mixes 
CDW, CF, A, and rj types of order. The state |C 2 ) (identified with the CDW' phase) mixes 
CDW, PI, and SF types of order. The state (L^) which is connected with the PS phase is 
also a mixture of dSC and rj pairing. Finally the \W) state mixes SDW, dSC, CF, and PI 
types of order in complete agreement with the results for the susceptibilities in the two-patch 
approach. 

Although we do not analyze in this section the contribution of the whole FS, one can 
expect that general features learned from the simple two-patch analysis hold in this case 
too. Namely, we expect that the eigenstates of the whole Hamiltonian do not coincide 
with the eigenstates of operators, corresponding to different order parameters. Although 
the eigenvectors of the effective Hamiltonian corresponding to the contribution of the whole 
FS have a much more complicated form, we expect that the combinations of the order 
parameters listed in Table 2 which belong to the same eigenstate remain unchanged. This 
can be also seen from the analysis of simultaneously diverging susceptibilities in the many- 
patch approach. At the same time, one can not safely argue whether or not the discussed 
states support the coexistence of long-range orders, or whether part of the abovementioned 
orders are only short-range. Resolving these possibilities requires a strong-coupling analysis 
of the problem which is beyond the validity of one-loop renormalizat ion-group approach. 
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V. SUMMARY AND CONCLUSIONS 



We considered the phase diagrams of the extended U-V-J Hubbard model within mean- 
field approximation and from RG approaches. The extended Hubbard model has a very 
complex phase diagram with various types of orders. The mean-field approximation pro- 
vides only a rough phase diagram of the system. At positive J it leads mainly to the 
coexistence of antiferromagnetism with a small amplitude of rf-wave superconductivity. At 
negative J mean-field theory predicts ferromagnetic order. For large enough positive V these 
phases are replaced by CDW order, while for large negative V the ground state is a rf-wave 
superconductor. Charge- and spin- flux phases, as well as their coexistence with rf-wave 
superconductivity are never energetically stable within the mean-field approximation. 

The RG analyses at van Hove band fillings leads to substantially reacher phase diagrams 
(Figs. 4-6). Instabilities towards SDW, CDW, dSC, CF, SF, F order are possible in different 
parameter regimes. Importantly, for J > and not too large \ V\ SDW, CF, dSC orders have 
comparable susceptibilities, which signals their close competition in this parameter region. 
At the same time, the tendency towards the formation of charge-flux is greatly suppressed 
in the many-patch approach in comparison with the two-patch results, and it becomes the 
leading instability only in the restricted parameter range. With increasing \t'\ the rf-wave 
superconducting state becomes more preferable. Tendencies towards both spin-independent 
and spin-dependent Pomeranchuk instabilities are outside of the weak-coupling region of the 
phase diagrams where the RG analysis is applicable. 

The symmetry analysis shows that most of the phases obtained within the RG approach 
are in fact a mixture of long- or short-range orders of different types. Resolving between 
short- and long-range orders for these phases requires a strong-coupling analysis of the 
problem which can not be performed within the RG approach. 

Comparing the obtained results with the phase diagram of the ID version of the U-V-J 
Hubbard model [18], we observe that the phase diagram of ID model is similar to the results 
of the two-patch approach at half-filling. Charge- and spin-flux phases replace the staggered 
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bond-order-wave phase of the ID model. The CDW, CDW phases appear in two dimensions 
in the same parameter region as the CDW phase in the ID case. For J > the SDW phase 
of the ID system is partially substituted by <i-wave superconductivity. At the same time, 
the results of two-patch approach are strongly changed at half filling by the contribution of 
the whole Fermi surface and become closer to the mean-field phase diagram, rather than to 
ID results. 

The physically most relevant regime of the 2D U-V-J Hubbard model is contained in the 
parameter range < V < U/A, < J < U/2. The corresponding areas are shaded in Figs. 
4-6. With decreasing van Hove filling we observe the following qualitative changes in the 
shaded area of the phase diagram: at half-filling SDW order dominates. For small doping 
5=1 — iivh = 0.08, we encounter a very complicated situation, in which SDW, charge-flux 
and d-wave superconductivity ordering tendencies are all simultaneously strong, although 
SDW order still dominates. In this regime the ground state structure is very sensitive to 
small parameter changes leading to a variety of possible phase transitions and the possibility 
for coexisting phases. Finally, at larger doping 5 = 0.28 the picture becomes simpler again 
and only <i-wave superconductivity becomes the leading instability. 

Under the abovementioned circumstances of strong competition of different order pa- 
rameters, the self-energy effects which are not accounted for in the one-loop RG approach 
can become crucial. Therefore, for a final conclusion about the possibility of a charge-flux 
phase in a broad parameter range of the one-band U-V-J model, the analysis of the two-loop 
contributions to the RG schemes remains to be performed. Another topic for future work 
will be to investigate to what extent the results for the one-band model U — V — J model in 
a specific parameter range allow implications for the competing orders in cuprate materials. 
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Tables 



k 


E k ({ gi }) 


State 


Corresp. phases 


1. 


-92 - g± 


\Vi) = 


T,T), \\ 


/ 2) = |T,I) + |!,T), |V3) = |U) 


F, SF 


2. 





\Ui) = 


o,o) + | 


U,U>; |^2> = |0,0>-|U,U> 


PS 


3. 


2ffi -92-94 


\Ci) = 


T, I) - 


I!,!) 


CDW, A 


4. 


9i ~ 2# 2 + 53 


\C 2 ) = 


U,o) + |o,u) 


CDW 


5. 


ffi - 2# 2 - 93 


\W) = 


-1 U,o) + |o,u) 


SDW, CF, dSC, PI 



Table 1. The eigenvalues E k ({gj}) and eigenstates \k) of the effective Hamiltonian, Eq. (36). 
The right column represents the correspondence with the phases in the phase diagrams (see 
text). 



State 


|SDW) 


|dSC) 


|CF) 


|SF) 


|CDW) 


|r> 


\A) 


\v) 


(v*\ 


l/y/2 








i/V2 








l/\/2 





(Ci| 













-1/V2 





-1/V2 


-1/V2 


(C 2 | 











1/V2 


1/v^ 


l/v/2 













-1/V2 

















-1/V2 


(W\ 


1/V2 


-1/^2 


1/v^ 








-1/V2 









Table 2. The scalar products C mfe = (k\m) between the eigenstates of the effective Hamil- 
tonian, Eq. (36) (see also Table 1) and the eigenstates of the operators in Eqs. (2) and 
(24). 
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Figure captions 



Fig. 1. The mean-field phase diagram for U = At, t'/t = 0.1 and van Hove band filling 
n = 0.92. SDW and CDW denote the spin- and charge- density wave phases, dSC the <i-wave 
superconducting phase, SDW+rfSC marks phase coexistence, and F is the ferromagnetic 
phase. Solid and dashed lines correspond to first and second order transitions, respectively. 

Fig. 2. The Fermi surface at van Hove band fillings: t' — and n — 1 (solid line), 
t'/t = 0.1 and n = 0.92 (long-dashed line), and t'/t = 0.3 and n = 0.72 (short-dashed line), 
A and B are van Hove points. 

Fig. 3. Diagrammatic representation for the many-patch RG equations, Eq. (28). Lines 
drawn through the vertices show the direction of spin conservation. Diagrams are drawn in 
the same order as the respective terms in Eq. (28). The cutting dash at the propagator lines 
means the derivative with respect to T (for brevity we indicate only the derivative of one of 
the propagators, the same diagrams with derivatives of another propagator are included as 
well). 

Fig. 4. Phase diagram as obtained from two- and many-patch RG analyses for U = 2t, 
t' = 0, and n — 1. Bold lines bound the weak-coupling region of the phase diagram, where 
the RG approach is applicable. Solid lines correspond to the phase boundaries obtained from 
the two-patch analysis. CF and SF denote the charge- and spin-flux phases, respectively, 
PS is the interaction-induced phase separated state. F, PI, and A denote the possibilities 
for ferromagnetism (F), spontaneous spin-independent (PI) and spin-dependent (A) defor- 
mations of the FS in the corresponding strong coupling regions. The other notations are 
the same as in Fig.l. The shaded region is the frustrated area with the critical scaling 
parameter of the two-patch approach A c > 20. The captions at the phase boundary lines 
denote the symmetry of the Hamiltonian (16) on these lines with respect to the rotation 
operators R v , R n , R T or Ra and the corresponding ordered states which become equivalent. 
The symbols correspond to the results of the many-patch RG approach: circles correspond 
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to SDW, diamonds to dSC, squares to SF, triangles to F, crosses to PS, and stars to CDW 
phase. The open symbols denote the frustrated regime in many-patch calculations when no 
divergence of the coupling constants was obtained for At = ln(t/T)/2 < 4. The type of the 
corresponding ordering tendency in this case was determined by the largest susceptibility at 
\ T = 4. 

Fig. 5. Phase diagram as obtained from the RG calculations for U = 2t, t'/t = 0.1, 
and n = 0.92. Solid and dashed lines correspond to the phase boundaries obtained from 
the two-patch analysis and separate the phases with respect to the scaling behavior of the 
coupling constants (see text). The half-filled circle corresponds a to charge-flux instability 
in the many-patch RG analysis. Thin crosses bound the bottom of the region where the 
charge-flux fluctuations become substantial, xcf > (2/3)xsdw- Other notations are the 
same as in Fig. 4. 

Fig. 6. Phase diagram as obtained from the RG calculations for U = 2t, t'/t = 0.3, and 
n = 0.72. Thin crosses bound the region where charge-flux fluctuations become substantial, 
Xcf > (2/3)max(xg DW , Xdsc)- Other notations are the same as in Figs. 4,5. 

Fig. 7. Scale dependence of the leading order-parameter inverse susceptibilities in the 
two-patch RG approach for U = 2t, t'/t = 0.1, and n = 0.92: a) V = J = 0; b) V = 0, 
J = 0.2U; c) V = 0.5U, J = 0. The vertical dot-dashed lines mark the scaling parameter A* 
where the largest absolute value of the coupling constants \gi\ is unity. 
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